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On Whitham's Averaging Method 

A. Ya. Maltsev and M. V. Pavlov* UDC 517.944 

Introduction. Whitham's averaged equations [1] for a nonlinear evolution system describe slow mod- 
ulations of parameters on a family of periodic traveling wave solutions (or on families of multiphase 
quasiperiodic solutions, which are so far knows to exist only for integrable equations) and are a system of 
\ hydrodynamic type [2, 3], that is, of the form 

! Vr = V}(U)U i x , i,j = l,...,N (1) 

(we consider only the spatially one-dimensional case), where U — (U 1 , . . . , U N ) . 

The original evolution system is usually Lagrangian or Hamiltonian, and this property is then inher- 
ited by the equations of slow modulations. A Hamiltonian theory of systems (1) was constructed by 
B. A. Dubrovin and S. P. Novikov [2] and, then, successfully used by S. P. Tsarev [7] to integrate Hamil- 
tonian systems reducible to the diagonal form 

W T = V i {U)U i x , i = l,..., N. (2) 

In the present paper we give an exposition and justification of some topics in Whitham's averaging 
theory. Namely, we consider classical single-phase Whitham averaging for a Lagrangian system in which 
some of the fields (denoted by u % ) have a direct physical meaning, whereas for the other fields (denoted 
by tp a ) only the ^-derivatives (x) are physically meaningful. The averaging is carried out on a family 
of periodic traveling waves of the system rewritten in the coordinates u l (x) , q a (x) = </?"(x) . To average 
such systems directly on the basis of the Lagrangian formalism, Whitham [1] proposed the pseudophase 
method, which permits one to obtain the averaged equations in the Lagrangian form and then to con- 
struct the corresponding Hamiltonian structure [6] (for the case in which pseudophases are lacking, the 
Hamiltonian formalism for the Whitham equations was constructed by Hayes [11]). However, there is an 
\Q ' alternative approach. It is based on the fact that in the variables (u l (x), q a (x)) the system is Hamiltonian 

and possesses conservation laws of the form SI = TC X , which correspond to the energy and momentum 
conservation laws and to the annihilators of the Poisson bracket. By using these conservation laws, we 
can rewrite the equations of slow modulations in the form 



in 



o 



(S i ) T = (R i ) X) (3) 

where (...) stands for the averaging on the family of traveling waves [1]. The Hamiltonian structure 
of Eqs. (3) can be obtained by averaging the original Hamiltonian structure by the Dubrovin-Novikov 



X 



method [2, 3], but the general proof of the Jacobi identity for the averaged Poisson bracket is lacking [4]. 
In the present paper we prove that 1) both methods lead to the same equations of slow modulations (this 
was proved in [4] by the WKB method) and 2) the Poisson bracket averaged by the Dubrovin-Novikov 
method coincides with that obtained from the averaged Lagrangian formalism (thus, the Jacobi identity 
for the Dubrovin-Novikov bracket is proved in this case). Furthermore, we prove that if the system has 
additional conservation laws of the form St = R x , then this averaging results in equations consistent 
with (3). 

1. Averaging of conservation laws and the pseudophase method. Let a q (x,t) be the field 
functions, and let the action 



S = J L[cP(x, t)] dt = JjL<a', of, <4, aj t , <,<&,•••) dxdt, 



(4) 
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q = 1 , . . . , N , be given. 

Suppose that an action S G T m : (a q ) — > (a 9 ) of an m-dimensional Abelian group T m (not necessarily 
compact) on the manifold of fields is given, the orbits of this action are m-dimensional, and the Lagrangian 
is invariant, that is, 



H(a q ,a q t ,a q ,...)=L(d q ,a q ,a q ,...). 



(5) 



By Noether's theorem, each symmetry of the Lagrangian (including translational invariance and inde- 
pendence of time) yields a conservation law. It is easy to verify that these laws have the form 



E\ q dL 



-n>l q 

(conservation of energy) and 



n>l 



+ [...]tt + [---]xx + [---]xt = 



(6) 



<9£ 



1)4 ' x da q nt 



E n E a! 



„ an 

nx dal x 



+ [...} tt + [...} xx + [...} xt = (7) 



(conservation of momentum); here a q nt = d n a q /dt n , = d n a q /dx n . 

Before writing out the other conservation laws, let us pass to new field variables, (a q ) — > (u l , <p a ) , 
where u l are constant on the orbits of T m and ip a are parameters on these orbits (and coincide with the 
parameters on the group itself). 

In these variables the other conservation laws acquire the form 



" <9£ " 




' da ' 




+ 

t 





[ ■ ■ ■ ]tt + [ • • • \xx + [ ■ • • ]a 



a = 1, . . . , m. 



The equations of motion have the form 

dL d dL d dL 
da q dt da q t dx da% 



+ ■ 



0. 



1,...,N. 



(8) 



(9) 



Assume that the Lagrangian is nondegenerate. Then we can proceed to the Hamiltonian formalism by 
setting (the method of Ostrogradskii applied to the case of fields) 



tfi(x) = if a (x), ... , q% a (x) = ^( na _i) t (x), 
r\(x)=u i (x), r hi {x) = u\ n _ 1)t {x) . 



(10) 
(11) 



(Here n a and are the highest orders of time derivatives of (p a and u % , respectively, occurring in the 
Lagrangian. Since we can perform linear transformations of the coordinates (p a and u 1 (separately), it 
follows that in the generic case all n a (and all n,) are the same.) Let us introduce the momenta 



P?(*) = 



5L 



S\(X) = 



6<p?(x) 
SL 



. . + (_l)«o 



5L 



(n a -l)t 



(12) 



8<Pn a t(x) ' 

SL 



*a 4 (») = 



8u\{x) 

5L 

Su lt( x ) ' 



-••• + (-1) 



71,-1 



SL 



5u kt( x ) 



(j»i-i)t 



(13) 



where 



SL 8L , . d dL , , 



The equations of motion (9) in the variables (x) , r l (x) , p"^ (x) , s l (x) are Hamiltonian. The 
Hamiltonian is equal to 



H = f \ E p" a (*) c (*) + E 4 (*) 4 (*) 



cfo-L, (14) 



and the Poisson bracket has the form 

(CW.^W) - <5 Q/3 <W(x - j,), {^(x), s^(y)} = ^S^S(x y) (15) 

(all other brackets are zero). 

More generally, if the conservation law corresponding to the symmetry of the Lagrangian with respect 
to a one-parameter transformation family can be put in the form P t = Q x (perhaps, ambiguously, because 
of the presence of mixed t, ^-derivatives, as is the case in (6), (7), and (8)), then the functional J P dx 
(defined unambiguously) generates this transformation family in the Hamiltonian structure (15). 

Remark. Very frequently, an originally degenerate Lagrangian becomes nondegenerate after a suitable 
rotation in the plane (x, t) (for example, this is the case for the Korteweg-de Vries (KdV) and nonlinear 
Schrodinger (NLS) equations). 

Let us describe the set of functions to be considered in the sequel. Let a(x) satisfy the conditions 

3fc>0, SeT m Vx a(x + 2n/k) = Sa(x), (16) 

that is, in the variables (u l , ip a ) we have 

u l (x + 2Tr/k) = u\x), <p a (x + 2ir/k) = p a (x)+fi a , (17) 

and let the initial data be given in the form 

Vx ul lt (x + 2n/k) = ul lt (x), ^ t (x + 2n/k) = ^ lt (x), n>l. (18) 

Since £ depends only on the derivatives of ip a , the evolution equations (9) preserve each of the families 
(17) and property (18). 

In terms of the Hamiltonian variables conditions (17) and (18) acquire the form 

q?(x + 27r/k) = q°(x)+ f i a , q«Jx + 2ir /k) = q«Jx) , v a >\, (19) 
rj li (x + 2n/k) = rl li (x), p« a (x + 2n / k) = p? a (x) , 4 (x + 2w/k) = 4 (*) . (20) 

On the family (19), (20) we seek extremals of all possible functionals of the form 

X H H + A P P + E A « /a , x h, Xp, X a = const, (21) 

where H is the Hamiltonian, P is the integral of momentum, and I a are the functionals that generate 
the shifts of the corresponding angles ip a . Thus, we consider the equation 



X H H + X P P + E AcJ" 



= . (22) 



We assume that for any k , \i a , Xh , Xp , and X a the solution to Eq. (22) is unique modulo a translation 
along the x-axis and the action of a transformation in T m . 
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Comment. If for any k > 0, fi a , Xh , Ap, and A a there exists a solution a 



(k/n,^n,X) 



(k.p.xjO*)' then 

(x) is also a solution for any positive integer n and to each tuple (k, p., x) there corresponds 
countably many such solutions. Accordingly, we assume that our family of solutions has just this very 
form, and to each a^ k _ ^ (x) we assign the greatest admissible k and the smallest admissible p, . 

The solutions to Eq. (22) depend on the parameters k, [i a , Xp/Xr , X a /Xn , , and #o > where y>g 
and 9 are the initial phases corresponding to the action of T m and to the translations along the x-axis. 
The averaging affects the 2m + 2 parameters k, fi a , Xp/Xh , and A a /Ajf ; the averaged variables do not 
depend on c^g and #o > so that these last variables do not occur in the equations of slow modulations. 

Let us transform system (14), (15) into a new Hamiltonian system by setting 



<f(x) = fi x (x), p a (x)=p1(x) 



(23) 



(this substitution is well-defined, since the densities of the functional H , P , and I a depend only on the 
derivatives of q" with respect to x). 
The new Poisson bracket has the form 



{< (*) , < (</)} - - y) , ^ , ^ > 2, 

{iW.if!/)} = <^'<W0r - y), {q a (x) lP ^y)} = S^S'(x - y) 



(24) 
(25) 



(all other brackets are zero). 

Thus, J q a {x) dx and J p /3 (x) dx are annihilators of the bracket (24), (25), so that every Hamiltonian 
system with translation-invariant Hamiltonian has 2m + 2 first integrals. 

Lemma 1. The family of periodic traveling waves of the Hamiltonian system (14), (23)-(25) can be 
obtained from the family (19), (20), (22) by factorization with respect to the initial phases </?q . 

PROOF. It is easy to see that conditions (19) and (20) are equivalent to the 27r//c-periodicity of the 
functions occurring in (23)-(25). 

It follows from conditions (22) that qf t = — (Ap/A#)gf x — A Q /A# and that £ t = — {Xp/Xh)£ x for the 
other variables. By differentiating the first equation with respect to x and by substituting q" x = q a (x) 
wc obtain qf — — (Ap/A#)<7"; thus, after factorization with respect to the initial phases tp^ we obtain 
exactly the family of periodic traveling waves of system (14), (23)-(25). The lemma is proved. 

COROLLARY. The family of periodic traveling waves of system (14), (23)— (25) depends on (2m + 2) 
parameters (not including the initial phase 9q). The averaging of the conservation laws for the energy, 
momentum, and 2m annihilators of the bracket (24), (25), yields the Whitham equations of slow modula- 
tions. 



Since 



dL 



SH SH _ dL 



(26) 



8tf y " 5q«(x) 

modulo total derivatives with respect to x and t and the energy and momentum conservation laws have 
the form (6), (7), we see that the equations of slow modulations have the form 

3£U d 



V 



n>l 



E 



i dL 
nt d^ 



71 L 



E 



<Pnt 



+ 



E 

n>l 



dL 

" ( "- 1)a; <* Out. 



EME 



n>l 



dL 



dL 



E 

i 

E ' 



+E 



dL 



0, (27) 



x 



n>l 



E 



dL 

duL 



+ 



Vr. 



dL 



(L) 



0, (28) 



x 



9 



Biff I rp 

(Or = 



6H 



8p«(x)/ x 



(29) 
(30) 



where (...) denotes the averaging on the family (19), (20), (22) (or, which is the same, on the family of 
periodic traveling waves of system (14), (23)-(25)). 

It is easy to see that Eqs. (29) correspond to the conservation laws (8). 

Let us proceed to the averaging method for the Lagrangian. 

Following Whitham, we seek solutions to Eqs. (9) in the form 



u l {x,t) = $ l (6>), <p a (x,t) = y a (9)+e c 



(31) 



where 9 = kx + uit is the phase, e a = f3 a x + j a t are pseudophases, and <!>*(#) and ^ a (9) are 27r-periodic 
functions. 
Next, we set 

£(«Si4,4,...,^^,^x*.--^ 

where = $*(#, k, u, [3 a , j a ) , * Q = ^ a (9, k,ui, f3 a , •y a ) , construct the averaged Lagrangian 



2* Jo 



£,(0, k, uj, P a ,r) d9 = C(k, co,(3 a ^ a ) = L(9 X ,9 T , s x ,s%), 



and obtain the averaged equations 

[Ze x ]x + [Zg T ] T = 0, [l £% ]x + [Z £ *] T + 0, 

or, in the variables k , uj , [3 a , -f a , 

[£w]t + [&k]x = 0, k T =aj x , 
[Z r .] T + [Z f i*]x = 0, Pt=1x- 



(32) 

(33) 

(34) 
(35) 



Since the Lagrangian £ is independent of time and translation- invariant, Eqs. (34), (35) admit energy 
and momentum conservation laws. On replacing Eqs. (34) by these laws, we obtain the averaged system 
in the form 



a 

kZ^ + ^^Z. 



a 

kZ + ^2 p a Zp~ - z 



0. 



T 

[V] T + [£/3«k = 0, 

Pt = 7x ■ 



(36) 

(37) 

(38) 
(39) 



Lemma 2. TTie solution family (31) coincides with the family (19), (20), (22), and moreover, [3 a = 
(O, 7 a --(Ap/A H )(9 a ) 1 A Q /A H . 

PROOF. The coincidence of these solution families is evident, and (22) implies that 

f3 a = (k/2n)[p a (x + 27r/fc) - <p<*(x)] = «) = (q") ; 
similarly, we have -f a — (ipf) — (qf t ) , but, by virtue of (22), (22) qf t = — (Ap/Ajj)gf x — X a /^H , that is, 
7 Q - (-(Ap/A ff )(fe - \ a /\ H ) = -(Xp/Xh) (q a ) - K/\h ■ 
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The lemma is proved. 

It is also obvious that the averaging with respect to is equivalent to the averaging with respect to x 
and t. 

Theorem 1. Equations (27)-(30) coincide with Eqs. (36)-(39). 
Proof. Let us introduce the functions 

Z-— £"- dL L n - dL L n -— L n - 9L r>>^ 

l ~du^ Li ~du^ L «~d^ x > ^~du^ La ~d^ ~ 
For all these functions we set 

LI = HI , w* j , k% , . . . , + 7 « , + (3 a , u:H a Q e ,k 2 ^ e ,...), 
where & = ¥(6, k, u},0,j), = * Q (6», k, us, j3, 7), and define the action 
~ 1 f 2jT 

S= 2^Jo ^ i ^^e,m,...,^ a e + 1 a ,k^ + p a ,w 2 ^ ,kH a ee ,...)de. 



(Obviously, 5= (XL).) 
It is easy to verify that 



6&(0) ' <J* tt (0) 
on the solution family (31). Furthermore, it is obvious that 



ss ss 

0, 7^T=0 (40) 



^ = E n ( E <" n_1 <*w> + E 

n>l ^ » a y 

n>l i a 



(here notation such as £ and (XL) is used for the averaged variables). 
By (40), we have 



n>l ^ j a ' 

E * = e »(e fcn_i <*n^r> + e * n-i »£s>) . 



n>l 
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It follows from Eqs. (27)-(30) and from Lemma 2 that 



n>l 



n>l ^ i ' n>l ck ' a 



E "( E ("(n-l)*,* + E (^n-1) 



n>l 



= E ^(E( fc " _1 ^ £ ") + E( fc " _1 ^ £ «)) + E(^«) 

n>l ^ i a 'a 

= uZk + ^^Zpc. ; 

a 

thus, (27) coincides with (36), and 



= E»(E< wn_lfc *w> + E< wn ~ lfc *^s> ) + E<^> 

n>l i a 

a 

e<e(<.^)+e(^))-<,> 

= E«fE( fcn<i> «^")) +E f E( fcn *n«) +E(/ ?Q ^)- w 



n>l 



n>l x i 7 n>l x a 



which means that (28) coincides with (37). 
Furthermore, 



that is, (29) coincides with (38). 
Now, by (22), 

mn) " ( " £ WW - P » - S & " ( " £ * " &) ~£ <«"> - & ' 

and (30) coincides with (39) by virtue of Lemma 2. 
The theorem is proved. 

2. Additional conservation laws. Let us now give an independent proof of the fact that the 
averaging of the conservation laws corresponding to energy, momentum, and I a on the family (19), (20), 
(22) and the equations (3% — -f x (which correspond to the averaging of the annihilators J q a {x) dx of the 
Poisson bracket (24), (25)) imply the "conservation of waves" kr = <^x ■ Furthermore, we shall prove that 
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if the Hamiltonian system (14), (15) has additional conservation laws St — R x corresponding to some 
fluxes on the space of a q {x) which preserve the Lagrangian and commute with the fluxes generated by the 
integral of momentum and by I a (in this case their densities and the fluxes themselves are independent of 
tp a and, in the variables used in (24), (25), correspond to conservation laws for system (14), (24), (25)), 
then the averaging of these additional conservation laws yields equations consistent with (27)-(30). 

Theorem 2. The averaging of the conservation laws for the energy, momentum, and 2m annihilators 
of the Hamiltonian system (14), (24), (25) on the family (19), (20), (22) implies the conservation of waves 
kp = ujx ■ 

Proof. In accordance with the preceding, we consider the equations 



X H H + X P P + J2 ^1° 



= 



on the phase space of functions that satisfy the conditions 

u l (x + 2ir/k) = u\x), <p a (x + 2ir/k) = tp a (x) + 2ir/3 a /k . 

We assume that the variations Su l (x) , 5ip a (x) are uniformly bounded, that is, u l (x) + Su l (x) , ip a (x) 
8tp a (x) belong to the same family (15) as u l (x), ip a (x). Thus, 



5ai(x 



Any functional that commutes with the Hamiltonian, with the momentum, and with I a leaves each of 
the families (19), (20) invariant. 

Let U = (U 1 , . . . , JJ 2m+2 ) be a collection of parameters (excluding the phases and 8q) on the 
solution family (19), (20), (22). Let £ be a tangent vector to the level surface k = const, (3 a — const in 
the space with coordinates (U 1 , . . . , JJ 2m+2 ) . Then the function variations correponding to the translation 
along £ are uniformly bounded. Set 



H = j'J>Hdx, P^Jypdx, I a = J'J> a dx 



(the corresponding conservation laws have the form = (Jh)x , i^p)t = (Jp)x , ('Pa)t — (J a )x)- By 

(22), we have 

\ H {u)d^ H ) + \p(V)d£(9 P ) + Y / X «(u)d i (y a ) = o 

a 

for any £ such that dgk — and d^[3 a = . This equation is equivalent to the relation 

X H (U)d(y H ) + Xp(U) d(V P ) + K(U)d(V a ) - n(U) dk + Y, V*(U) d(3 a 

a a 

for some functions (i(U) and fj, a (U) . 

As was already shown, under conditions (19) and (20) the solutions to (22) have the form (31), that is, 

u l {x,t) = <5> l {0), (f a (x,t) = ^ a (9)+e a , 6 = kx + ojt, e a = [3 a x + j a t, 

where & l (9), \&"(0) are 27r-periodic functions. 

In the original Lagrangian formalism let us make the rotation by an angle \ m the plane (x,t): 

x = xcosx — t sinx, t — x sin \ + t cos \ . 

The passage to the Hamiltonian structure will now give a new Poisson bracket, in which the translations 
along the "old" t-axis will be generated by the functional J (Vh cos x~ Jh sin \) dx , the translations along 
the old coordinate x by the functional /(3\p cos x — Jp sin x) dx , and the action of T m by the functional 
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/ (lP a cos x ~ J a sin x) dx , where all densities and fluxes of the conservation laws are expressed via the new 
coordinates a 9 (f) . (Indeed, 7 t = J x -<=^> (J'cosx — Jsinx) t " = (IPsinx + Jcosx)^ .) 

However, the solution family (31) remains the same, and moreover, k — fccosx — wsinx, $ a = 
(3 a cos x — 7" sin X ; w = fc sin x + ^ cos x , 7 a = /3 Q sin x + 7° cos x • 

By writing out the conditions 



AS 



X H J y H dx + \ P j ? P dx + J2 A a J 



in the new variables for some function in the family (31), we obtain 



A( X )S 



cos x — Jh sin x) dx + Xp / (!Pp cos x — Jp sin x) da; 



= o, 



where A and A(x) are the old and the new Hamiltonian operators, respectively. 
Since A( X ) is nondegenerate, we obtain 

h J^ H cos X — Jh sinx) dx + Xp J^p cos x — Jp sin x) dx 

cos X- Ja sin x) (if 



0. 



Since the averaging on the family (31) in any direction in the plane (x,t) (except for 9 = const, 
e a = const) gives the same result, we can conclude, as above, that 

X H (U) {d ( 7 H ) cos x - d ( J H ) sin x ) + Ap (U) (d (5> P ) cos X - d ( J P ) sin X ) 

+ K(U)(d(J> a ) cosx - d(J a ) sinx) 



that is, 



= f i( X ,U)dk + J2^(x,U)d/3 a , 

a 

X H (U) (d C? H ) cos x - d (J H ) sin X ) + Ap (tf) (d (?p) cos x - d ( Jp> sin X ) 
+ ^ A Q (J?) (d (0> a ) cos x - d ( J Q ) sin x) 

= KX,U) (dk cos x - dui sin x) + (x > ^) (d/3 Q cos X ~ dj a sin x) • 

a 

For x = and x = 7r /2 we obtain, respectively, 

x H (U)d(y H ) + x P (U)d(y P ) + ]T xju) d(v a ) = M (o, to dfc + ]T m«(o, i7) , 

ck a 

A ff ([/) d(T ff ) + Xp(U)d(J P ) + £ X a (U)d(J a ) = h(tt/2, U) dco + ^ Mtt/2, tf) ^7°- 



(41) 

(42) 
(43) 



Let us multiply (42) by cosx and (43) by sinx and subtract the results from (41). We obtain 
[M(X, U) - m(0, t 7 )] cosxdfc - [/x(x, 17) - M 71 "/ 2 , U)\ sm X duj 

+ £>a(x, C7) - M(0, C7)] cosxd/3« - J>«(x, tf) - M«(V2, tf)] sinxd 7 « = . 
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Assuming that the differentials dk, dcu , df3 a , and dj a are linearly independent (on the space (U)), we ob- 
tain n(x, U) =ju(0, U) = /u(tt/2, U) = /jl(U), /j, a (x, U) = /z Q (0, U) = fi a {n/2, U) = fj,(U). Consequently, 
on the space (U) we have 

X H (U) d(? H ) + \p(U) d{f P ) + M^) d{y a ) = v{U) dk + J2^{u) dP a , 

a a 

X H (U)d(J H ) + Xp(U)d(J P ) + K(U)d{J a ) = fi(U) duj + J2 Ha(U) dj a ; 

these relations readily imply the statement of the theorem. The theorem is proved. 

Theorem 3. Suppose that for any tuple Xh , Xp , X a , k, (5 a system (22) has a unique solution in the 
class (19), (20) modulo translations along the x-axis and the action of I a . Let the Hamiltonian system 
(14), (15) have an additional first integral St — R x , corresponding to some flux on the space of a q (x) , 
and suppose that this flux commutes with the fluxes generated by H , P , and I a and preserves the action. 
Then Eqs. (27)-(30) imply that (S) T = (R)x ■ 

Proof. The flux generated by the functional J S dx preserves the solution family (19), (20), (22) and 
can at most generate a linear dependence of the phases tp^ , 6q on time on this family. Since the Poisson 
bracket (15) is nondegenerate, it follows that 



J Sdx + ({U) J 7>pdx + J2(»(U) J y a dx 



on the family (19), (20), (22) for some function ((U) and Ca(U) . Literally repeating the argument in the 
proof of Theorem 2, we see that 

d(S) + C(U)d{9 P ) + ]T Ca(U) d(9 a ) = V (U) dk + Y, Va(U) df3 a , 

a a 

d (R) + C(U) d (J P ) +J2UU)d(J a )=r)(U)du; + Y Va (U) d 7 a 

a a 

on the manifold (U) for some functions r](U) , rj a (U) . The desired assertion now follows readily from 
Theorem 2. The theorem is proved. 

3. Hamiltonian formalism. Equations (34), (35) are Hamiltonian with respect to the Poisson bracket 

{[3 a {X),Z^{Y)} = 5 a H'{X-Y), {k{X),Z UJ {Y)}=5'{X-Y) (44) 

(all other brackets are zero); the Hamiltonian is equal to 

H = J ^ W + ^ 7 Q £ 7 . -Z^jdX, (45) 

the integral of momentum 

P = J (kZ^ + Y^&j^dX 

generates the translation along the X-axis, and the functional / P a {X) dX , J £ 7 a (X) dX , / k(X) dX , 
/ &uj(X) dX are annihilators of the bracket (44) (see [6]). 

Here we shall prove that the bracket (44) coincides with the averaged Dubrovin-Novikov bracket [2, 3] 
for Eqs. (27)-(30). 

Let us describe the construction of the Dubrovin-Novikov bracket. Let the original evolution system 
a t (x)=K(a,a x ,...) = {a(x),H[a}}, a={a q ), (46) 
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be Hamiltonian with respect to the local translation-invariant field-theoretic Poisson bracket 

M 

{a«{x),a?{y)} = ^Bf (a(x), a x (x), . . .)5^(x - y), (47) 

fc=0 

and let the Hamiltonian H[a] be a local field functional, 

H[a] = J h(a(x), a x (x), . . . ) dx . (48) 
Suppose that system (46) has N pairwise commuting local integrals 



>] = Jr{a{x),a x {x),...)dx, i = l,...,N, {P,P} = 0. 



(49) 



(In our case, N — 2m + 2 and the integrals have the form J q a (x)dx, jp a (x)dx, jVpdx, and / Vndx.) 
Consider the pairwise brackets of the densities of the integrals (49) 

{T(a(x),a x (x),...),T(a{y),a y {y),...)} 

= J2A^(a(x),a x (x),...)5^(x-y), i, j = 1, . . . , 2m + 2 . (50) 

k>0 

By (49), Aq — d x Q % i . The averaged Dubrovin-Novikov bracket has the form (we denote (V 1 ) = U l ) 

{U i (X),W(Y)} = (A[ j )(U(X))S'(X — Y)+ diQ ^ U) U k x S(X - Y) 

= g^{U{X))5'{X -Y) + tfi{U{X))U k x 5{X Y), (51) 

where (...) stands for the averaging on the family of periodic traveling waves of the Hamiltonian system 
(46), determined by the conditions 



k J "Ph dx — uj J Tpdx + ^ (^i a j q a (x) dx + X a J p a (x) dx^j 



0, (52) 



where k is the wave number, u is the frequency, and [i a and A Q are arbitrary constants. Just as in the 
proof of Theorem 2, it follows that 

k(U)d{T H ) - u;(U)d(yp) + J2(^(U)d(q a ) + X a (U)d(p a )) = M (C7) dk (53) 

a 

for some function n(U) . 

The theory of the brackets (51) is closely related to Riemannian geometry; in particular, their skew 
symmetry implies that 

g ij =g ji , q+^ = dg^/dU\ (54) 

and it follows from the Leibniz identity that under invertible smooth changes U l — > U l (U) of the field 
variables the functions g tj behave as the contravariant components of a metric tensor, and the functions 
Tj k = —gj S bf behave as the components of the differential-geometric connection consistent (by (54)) with 
the metric. 

If g 1 ^ is nondegenerate, then the Jacobi identity for the bracket (51) is equivalent to the symmetry of 
the connection r* fe and to the vanishing of the curvature tensor of the metric, Rj^i = (see [2, 3]). 

However, the survey [3] does not contain the proof of the Jacobi identity for the bracket (51) obtained 
by averaging (see [4]). In other words, it is not proved that the connection T l - k is symmetric and that the 

curvature tensor of the metric g lJ (U) is zero. By proving that the Dubrovin-Novikov bracket coincides 
with the bracket (44) rewritten in the variables ((q a ) , (p a ) , (T'p), (^h)) (as was shown in the proof of 
Theorem 1, we have (q a ) = (3 a and (p a ) — £ 7 a), we at the same time prove that for the case in question 
the Dubrovin-Novikov bracket satisfies the Jacobi identity. 
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Theorem 4. The Dubrovin-Novikov bracket (51) coincides with the bracket (44) rewritten in the vari- 
ables ((q a ),( P a ),(y P ),(y H )). 

Proof. Lot the bracket (44) in the variables ((q a ) , (p a ) , (3> P ), (Pjj)) have the form 

{U\X),W(Y)} = g^(U(X))S'(X — Y) + ~b^(U(X))U x S(X - Y) . (55) 

It is easy to verify [3] that the Dubrovin-Novikov bracket can be obtained as follows. With the bracket 
(47) we associate a bracket on the space of fields a q (x, X) by setting 

M 

{a«(x, X),aP(y, Y)} = (a(x, X),a x (x, X),.. .)S^(x - y)S(X - Y), (56) 

and then in the densities of the integrals - ^ (X) = J (a(x , X) , a x (x , X) , . . . ) dx and in the bracket (56) 
we replace the operator d x by d x + ed x , where e <C 1 is the ratio of the "rapid" scale to the "slow" scale 
(for example, the ratio of the wavelength to the typical length of variation of the parameters U). Then 
we evaluate the Poisson bracket of the resultant integrals P e {X) in the new Hamiltonian structure: 

{a?(x,X),ar(y,Y)} e 

M f 8 r) \ k 

= Y,Bf {a{x,X),a x {x, X)+ea x {x,X), ...){ — +e—\ 5(x - y)S(X — Y) . (57) 

k=o ^ ' 

Set 

h${X, Y) = {P £ (X),P £ (Y)} e = sh\ j (X, Y) + 0(e 2 ) 

and consider the restriction of the function h % l (X , Y) = htf{[a], X,Y) to the submanifold of functions 
a(x, X) such that for any fixed X these functions, regarded as a functions of x, belong to the family of 
traveling wave solutions to system (46). The function (3 3 *)(JC), 8 n (X) can serve as coordinates on this 
manifold, and the Dubrovin-Novikov bracket has the form 

{(V i ){X),(&{Y)} = h?{X,Y), 

where (X , Y) are the functions (X , Y) restricted to the cited submanifold (it is easy to see that 
they are independent of 9 (X)). It readily follows that in this Hamiltonian structure the functionals 



J(? H )dX, J(y P )dX, J(q a )dX, J(p a )dX 



(58) 



generate the averaged equations (27)-(30), the translation along the X-axis, and the "zero" fluxes, re- 
spectively on the space of fields U l (X) . 

As was shown in the proof of Theorem 1, (Vh) , (3V) , and (q a ) , (p a ) coincide, respectively, with 
the densities of the Hamiltonian, of the momentum, and of the annihilators of the Hamiltonian system 
(44), (45), and the fluxes generated by the functionals (58) by virtue of the Dubrovin-Novikov bracket 
coincide with the corresponding fluxes for the bracket (44) . Since the related equations in the coordinates 
(T){X) = U^X) have the form V\. = b ; k j (U)U x and XJ\. = b%(U)U x , respectively, it readily follows 

that b % £ = b % l and consequently, by (54), g tJ = g' L] + f lj , where f tJ = const. 
Furthermore, let 

r e (X) = r (X) + eIl(X) + O(e 2 ), 
{ ! { !, • "{ !: • ^: 2 i: 

then 

h\ j (X, Y) = {P (X),li(Y)} + {li(X),4(Y)} + {P (X),P (Y)} U 
and, by (52) (we set P = / 2m + 1 ; H = I 2m + 2 ) on the cited submanifold we have 

^r J (X)r J (r)^'(X,y)^0, 
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where (n) = (fi\ , . . . , /j, m , Ai, . . . , A TO , — w, k) . Since h % i(X,Y) has the form (51), we see that this 
relation is equivalent to the relation 



thus, (Ti(U)) = (fj,i(U), . . . , n m (U), \i(U), . . . , X m (U), —cu(U), k{U)) is an isotropic covector with re- 
spect to the metric g z; >(U) . In view of (53), we see that in the variables ((3 a , £ 7 q , k, C u ) this covector 
has the coordinates (0, . . . , 0, 0, . . . , 0, fi, 0) and is isotropic with respect to the metric (44). It follows 
that in the variables (U % ) the covector (ni(U), . . . , fj, m (U), Xi(U) , . . . , X m (U) , —u)(U), k(U)) is isotropic 
with respect to both the metric g 1 ^ (U) and the metric g l ^{U) , that is, it is isotropic with respect to pi . 
However, = const , and ^ a , \ a , u> , k can be arbitrary by assumption. Thus, any covector is isotropic 
with respect to / u ; that is, = 0. The theorem is proved. 

Example. The multicomponcnt NLS-type equation 



«4 + < x + n/*K = o, h = Y J \A\ i = i 

i 

after the substitution u % — Vufcx_p(i j v l dx) takes the form 



v, = 



w\ = -2(w i v i ) x . 
The system has the Hamiltonian structure 

{v i (x),w k (y)} = S ik S'(x-y), 

i \2 



(59) 
(60) 



*=/[e(-W-^ 2 M^ 



'IV 



dx, 



and the integral of momentum has the form P = J J2 i v l w l dx. 

The bracket has the form described in the preceding, and we can proceed to a nondegenerate Lagrangian 
formalism in the system after the substitution w l = q l x with allowance for the fact that v % = —q l t /{2q l x ) 
(thus, in accordance with the general procedure, we regard v l as generalized momenta and express them 
via q\). The system takes the Lagrangian form with the Lagrangian 



£ = E 



<l:r. 



Ml Ml 

We find the family of traveling waves of system (59), (60) from the conditions 





+ 


7< 2 y 


_2w i _ 


XX 


\2w l J 



We 



w 



cw x = -2(wV) a 



It follows that 





+ 




_2w\ 


X 


\2w z J 



cv 



that is, v 1 = (B l - cw 1 )/(2w l ) and 

r.2 



{v l ) 2 + V (J2 wk )+ At > cw t + 2w t v I ^B l 



+ 



2w l 



2w l 



\2w i 



We 



w 



0. 



(61) 



18 



System (61) has the first integral 

and is Hamiltonian with generalized momenta p l = w l x /(2w l ), Hamiltonian 

h = E (t - 4 + - V) 2 + iE? + Al A + v ( E - l ) 

and bracket {w 4 , p J } = <5 y . 

By applying the canonical transformation q' = 2Vvf , J? = Vw*p\ we obtain 

* . ? ( w + (ffi ) 2 + + Ml) + „( ? <£) , . 

Thus, the periodic traveling waves of system (59), (60) correspond to closed trajectories of a particle that 
moves in the n-dimensional space with the potential 

»»-?((SNfN^)V<?^) 

with all possible c, A 1 , B l . Assuming that (as usual) closed trajectories are isolated on each energy level, 
we obtain a (2n + 2)-parameter (excluding the initial phase shift 9q) family of traveling waves (or several 
such families) with the parameters c , A 1 , B % , E , on which we can average the 2n + 2 first integrals P , 
H , / v l dx, and / w l dx; all statements proved above are valid for this family. 

In closing, the authors express their gratitude to S. P. Novikov for his attention to the work, and also 
to I. M. Krichcvcr, E. V. Ferapontov, and O. I. Mokhov for useful remarks and fruitful discussion in the 
course of the present study. 

References 

1. G. Whitham, Linear and Nonlinear Waves, Wiley, New York (1974). 

2. B. A. Dubrovin and S. P. Novikov, "Hamiltonian formalism of one-dimensional systems of hydrody- 
namic type and the Bogolyubov- Whitham averaging method," Dokl. Akad. Nauk SSSR, 270, No. 4, 
781-785 (1983). 

3. B. A. Dubrovin and S. P. Novikov, "Hydrodynamics of weakly deformed soliton lattices. Differential 
geometry and Hamiltonian theory," Uspekhi Mat. Nauk, 44, No. 6 (270), 29-98 (1989). 

4. S. P. Novikov and A. Ya. Maltsev, "The Liouville form of averaged Poisson brackets," Uspekhi Mat. 
Nauk, 48, No. 1 (289), 155-156 (1993). 

5. M. V. Pavlov, "The Hamiltonian structure of Whitham's equations," Uspekhi Mat. Nauk, 49, No. 1 
(295), 219-220 (1994). 

6. M. V. Pavlov, "Double Lagrangian representation of the KdV equation," Dokl. Ross. Akad. Nauk, 
339, No. 2, 157-161 (1994). 

7. S. P. Tsarev, "On Poisson brackets and one-dimensional Hamiltonian systems of hydrodynamic type," 
Dokl. Akad. Nauk SSSR, 282, No. 3, 534-537 (1985). 

8. I. M. Krichcver, "Averaging method for two-dimensional 'integrable' equations," Funkts. Anal. Prilo- 
zhen., 22, No. 3, 37-52 (1988). 

9. I. M. Krichever, "Spectral theory of two-dimensional operators and its applications," Uspekhi Mat. 
Nauk, 44, No. 2, 121-184 (1989). 

10. M. J. Ablowitz and D. J. Benney, "The evolution of multiphase modes for nonlinear dispersive waves," 
Stud. Appl. Math., 49, No. 3, 225-238 (1970). 

11. W. D. Hayes, "Group velocity and nonlinear dispersive wave propagation," Proc. Roy. Soc. London, 
332, 199-221 (1973). 



19 



12. H. Flaschka, M. G. Forest, and D. W. McLaughlin, "Multiphase averaging and the inverse spectral 
solution of the Korteweg-de Vries equation," Comm. Pure Appl. Math., 33, No. 6, 739-784 (1980). 

13. J. C. Luke, "A perturbation method for nonlinear dispersive wave problems," Proc. Roy. Soc. London 
Scr. A, 292, No. 1430, 403-412 (1966). 

14. S. Yu. Dobrokhotov and V. P. Maslov, Finite-Gap Almost Periodic Solutions in the WKB Approxima- 
tion. Contemporary Problems in Mathematics [in Russian], Vol. 15, Itogi Nauki i Tekhniki, VINITI, 
Moscow (1980). 

Translated by V. E. Nazaikinskii 



20 



